Evolution of the accretion flow in Nova Muscae 1991 
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ABSTRACT 

We identify the X-ray reflected component in the Ginga spectra of Nova Muscae 
1991, a Black Hole transient system used as the prototype for the recent model of Esin, 
McClintock & Narayan (1997) based on advection dominated disk solutions. We see 
that the reflected spectrum is generally significantly relativistically smeared, and use 
this, together with the amplitude of reflection, to track the innermost extent of the 
accretion disk. The optically thick disk switches from being highly ionized to nearly 
neutral during the transition from high to low state, and the inner radius of the disk 
moves outwards during the low state decline. 

Qualitatively, this overall trend is compatible with Esin et al.'s model, but quantita- 
tively, the retreat of the inner disk during the high to low state transition is much slower 
than predicted. The hard (low state) spectra are not produced solely by an optically 
thin accretion flow: optically thick material within ~ 20 — 100 i? g is generally present. 

Subject headings: accretion, accretion disc - black hole physics - binaries: close - stars: 
individual (GS 1124-68) - X-ray: stars 

1. Introduction 

Black hole binary systems give one of the most direct ways in which to study the physics of 
accretion disks. There is no surface boundary layer or strong central magnetic field to disrupt the 
flow, and the orbital parameters are often well studied so that the inclination, mass and distance of 
the system are tightly constrained. Additionally, many of these systems (the Soft X-ray Transients, 
hereafter SXT) show dramatic outbursts where the luminosity rises rapidly from a very faint qui- 
escent state to one which is close to the Eddington limit, and then declines again over a period of 
months, giving a clear sequence of spectra as a function of mass accretion rate. The usual pattern 
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for such objects is for the outburst to be dominated by a soft component of temperature ~ 1 keV, 
with or without a (rather steep) power law tail, while the later stages of the decline show much 
harder power law spectra, extending out to 100-200 keV (see e.g. Tanaka & Shibazaki 1996). 

The "standard" accretion disk model developed by Shakura & Sunyaev (1973; hereafter SS) 
derives the accretion flow structure in the limit when the gravitational energy released is radiated 
locally in a geometrically thin disk. Such models give temperatures of order 1 keV for SXT in 
outburst, but are unable to explain the presence of a (hard or soft) power law tail to high energies. 
Either there are parts of the disk flow in a different configuration to that of SS, or there is some 
non-disk structure such as a corona powered by magnetic reconnection (e.g. Haardt, Maraschi & 
Ghisellini 1994). 

Recently, another solution of the accretion flow was postulated to explain the hard X-ray data. 
Below a critical accretion rate, m C rit) a stable, hot, optically thin, geometrically thick solution can 
be found if radial energy transport (advection) is included (see e.g. the review by Narayan 1997 
and the discussion by Chen et al. 1995 of how all the accretion flow solutions fit together). In the 
absence of an optically thick disk the hot electrons can only cool via cyclotron/synchrotron emission 
(on an internally generated magnetic field) , bremsstrahlung, or Compton scattering of the resultant 
spectra of these two processes. The lack of strong flux of soft photons in these models means that 
the X-ray spectra are typically hard. Such flows were proposed to explain the hard and very faint 
X-ray spectra seen from black holes candidates (BHC) in quiescence (Narayan, McClintock & Yi 
1996). These models were then extended by Esin, McClintock & Narayan (1997; hereafter EMN) 
to cover the whole range of luminosity seen in SXT. EMN assume that accretion takes place via an 
SS disk for the incoming material beyond some (large) truncation radius of r tr ~ 2 x 10 4 i? g (where 
R g = GM/c 2 ), but that interior to this the flow always jumps to the advective solution if the latter 
exists. During the outburst the mass accretion rate goes above that at which an advective flow 
can be sustained, so they assume that the disk adopts the SS solution, forming a cool flow down 
to Rm S = 6 R g and hence producing soft spectra. During the decline the mass accretion rate falls 
below m cr it, producing an intermediate stage as the SS disk retreats. The transition radius between 
the advective and SS flow progressively increases, and the low state proper is defined to be where 
the SS disk is again restricted to beyond 2 x 10 4 R g , although the spectrum is indistinguishable from 
this for transition radii > 200 R g . This change from an SS to advective disk solution takes place 
over a very small range in total luminosity, and is postulated to explain the soft-hard transition 
seen in the SXT. EMN tested their model on Ginga data of Nova Muscae 1991, and found that it 
could indeed reproduce the overall spectral evolution of the source. 

This model predicts large changes in the accretion geometry as a function of mass accretion rate 
(see Figure 1 of EMN) which can be observationally tested by spectral studies of X-ray reprocessing. 
X-rays interact with optically thick material to produce a reflected continuum due to Compton 
scattering and the iron Ka line due to fluorescence/recombination. Their normalizations determine 
the solid angle of the reprocessor as seen from the X-ray source while the shape of the reflected 
continuum below 10 keV, the line energy and its equivalent width are functions of ionization state 
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of the reprocessor. If the reflecting material is indeed a disk extending close to the black hole then 
the reprocessed spectrum is strongly smeared by the combination of high orbital velocities and 
strong gravity (Fabian et al. 1989; Laor 1991; Ross, Fabian & Brandt 1996). In the EMN model 
it is clear that as the source changes from high to low state then the amount of reflection should 
decrease, together with the amount of relativistic smearing, as the optically thick reflecting disk 
moves outwards. 

In this Letter we present results of spectral analysis of Ginga data of soft X-ray transient 
GS 1124-68 (Nova Muscae 1991) with the emphasis on detailed modeling of the effects of X- 
ray reprocessing. Our results confirm the prediction of EMN that the covering fraction of the 
reprocessing matter decreases as the source luminosity declines. However, the normalization of the 
reflected component and the amount of relativistic smearing are both significantly larger than the 
model predictions, indicating that the SS disk retreats outwards more slowly than predicted during 
the intermediate and low states. The hard (low state) spectra seen in the persistent BHC and 
SXT during decline are incompatible with the EMN scenario of an optically thin accretion flow 
extending out to ~ 10 4 i? g , and instead require that there is a transition from an optically thin to 
optically thick structure at ~ 20 — 100 R g . 

2. Data 

Observations of Nova Muscae 1991 and various aspects of data modeling were described in 
detail by Ebisawa et al. (1994). We re-extracted the data from the Ginga archive in Leicester 
and subtracted the background using a modified version (Zycki, Done h Smith 1997) of standard 
procedures (Hayashida et al. 1989). For our analysis we have selected six data sets: January 11th, 
May 18th, June 13th, June 21st, July 23rd, and September 3rd. These are the only datasets in 
which the hard tail is sufficiently developed to search for reflection, and where the instrument is 
accurately pointed (within 0.3° — 0.4°) at the source. 

3. Model for the hard X ray emission 

Our basic model of the X-ray spectrum is a sum of a primary power law and the reprocessed 
component which is a sum of the Compton-reflected continuum and the iron Ka line near 6.5 keV. 
The reflected continuum is computed using the angle-dependent Green functions of Magdziarz &: 
Zdziarski (1995), with the possibility of partial ionization (parameterised by ionization parameter, 
£ = Lx/nr 2 ) as in Done et al. (1992), all as implemented in the "pexriv" model in XSPEC version 
9.01 (Arnaud 1996). This is a very basic ionization code, balancing the photo-ionization from the 
power law irradiation with radiative recombination at a given fixed temperature, which we set to 
10 6 K. The self-consistent iron Ka line for this ionization state, computed using the Monte Carlo 
code described in Zycki & Czerny (1994), is then added to the continuum. 
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The relativistic effects expected from an accretion disk around a (non-rotating) black hole 
are simulated by convolving this total (line plus continuum) spectrum with the XSPEC "diskline" 
model. This is parameterised by the inner and outer radius of the disk, i?j n and R Q uti its inclination 
and the irradiation emissivity, F lir (r). The latter three are fixed at 10 4 i2 g , 60° (Orosz et al. 1996) 
and oc r -3 , respectively, while we fit for i?j n , the ionization state £, and the amplitude of the 
reprocessed component expressed as a solid angle subtended by the reprocessor from the X-ray 
source, £1, normalized to 2tt, i.e. / = Q/2n (see also Zycki et al. 1997). 

Our code also allows the ionization to vary as a function of radius, £(r) oc r' 3 , by splitting 
the disk into radial rings of constant ionization state. The relative covering fraction of each ring is 
determined by the emissivity and the relativistic corrections are given by the extent of each ring. 

4. Model fitting 

During its decline phase Nova Muscae 1991 went through a sequence of distinct "states" 
(Ebisawa et al. 1994; EMN). During and soon after the outburst the source was in the Very High 
State (VHS; January 11th data), with both a strong soft component and a power law tail. Two and 
a half months later that changed into High State (HS), where the soft component dominates, and 
further one and a half months later, in mid-May, the Intermediate State (IS) began in which the 
soft component dramatically decreases while the power law tail hardens and increases in intensity. 
This is then followed two months later by the Low State (LS; July 23rd and September 3rd data 
sets), where the hard power law dominates the luminosity output. An important point for the 
January (VHS) and May (HS/IS) data is the description of the soft thermal component since it 
dominates the spectrum below ~ 4keV. We model this by optically thick comptonization of soft 
seed photons (Titarchuk 1994; as implemented in the "comptt" model in XSPEC), since simple 
blackbody and disk blackbody models (with or without relativistic corrections) cannot give an 
adequate description of the observed spectral shape (see also Gierlihski et al. 1997b). The results 
for the June and July data are not sensitive to a choice of the soft component so here we use the 
multi-temperature disk blackbody model (Mitsuda et al. 1984). 

First, we have checked that the spectral features seen in the data at 5-8 keV, which we will 
attribute to X-ray reprocessing, cannot be accounted for by the effects of Comptonization, even 
if a broad distribution of soft seed photons is assumed. To this end we have modelled the May 
18th data using "comptt" for the soft component and we convolved it with the relativistic Green 
functions of Titarchuk (1994) to obtain the model for the hard component. Such a two component 
model fails to account for the data (x 2 = 656/25 dof ) and the residuals clearly indicate the presence 
of the reprocessed component. The situation is not improved if the soft component is modeled using 
the Green functions of Sunyaev & Titarchuk (1980) convolved with a blackbody. 

We first model the reprocessed component assuming its uniform ionization (£(r) = const; 
Table g). Reflect ion is significantly detected in all but the September 3rd dataset, where the 
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signal to noise is insufficient to obtain interesting constraints. Even in the VHS, where the soft 
component is very strong, reflection is significantly detected and its amount is quite well determined. 
It is strongly ionized and relativistically smeared (see Gierlihski et al. 1997b for a similar result for 
the soft state of Cyg X-l). At the beginning of the IS (May 18th data set) the reflector is still 
highly ionized, but as the importance of the soft component declines, the ionization state also drops 
(to £ ~ 10 where iron is predominantly ionized only to Fe VII-IX), and the power law hardens 
by Ar ~ 0.4 without a significant change in its total flux. Relativistic smearing of the reflection 
spectrum is significant in all the datasets except for July 23rd and September 3rd. The best-fit 
spectra for the May 18th and June 13th data sets are shown in Figure |l|. 

Figure || shows the derived confidence contours for R m and / during the IS and LS. Clearly the 
data require that the solid angle decreases as the source declines, and while there is no statistical 
requirement that the relativistic smearing decreases, the data are consistent with a monotonic 
increase in i?i n . 

To see how these results depend on any ionization structure, we repeated the fits for £(r) oc r~ L5 
and r -4,5 , as expected from a gas and radiation pressure dominated SS disk, respectively (with the 
assumed F^j- oc r~ 3 ). Only the contours for the May 18th dataset change, since this is the only 
one shown which has strong ionization. We include these results in Figure |^ by the dashed and 
dotted lines for the gas and radiation pressure dominated cases, respectively. Both allow much 
larger reflected fractions, since there is now a contribution from very highly ionized material at the 
inner disk edge which gives negligible spectral features. While this material is also strongly affected 
by the relativistic smearing, the lack of features means that there are no observational signatures 
until further out in the disk where the ionization drops, so allowing smaller values of i?i n . 

We have also checked that a different model for the soft component in the May data does 
not change the results concering the reprocessed component. Using the solution of Sunyaev & 
Titarchuk (1980) convolved with a blackbody spectrum to model the soft component we re-derived 
the confidence contours in the / vs. i?; n plane. The best fit and minimum values of / and i?i n agree 
to within 10% with previous results (solid curve in Fig. |^) while the maximum values are larger by 
20% when the Sunyaev &: Titarchuk solution is used. 

5. Discussion 

We can use the results of our spectral fitting to derive the geometry, and then compare this 
with that proposed by EMN in their scenario based on advective disk. Our derived limits on the 
inner radius of the accretion disk from the amount of relativistic smearing are dependent on the 
illumination emissivity, for which we assumed Fi„ oc r~ 3 . This should be a good approximation for 
a central quasi-spherical source even if there is some overlap between the cold disk and hot source 
(Poutanen, Krolik & Ryde 1997). Thus it seems likely that the radii derived in the previous section 
are reasonable unless the disk flares substantially, and/or there are significant non-local effects in 
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the energy generation and/or relativistic corrections are important. 

In the VHS, HS and beginning of the IS, where the spectrum contains a strong soft component, 
the EMN model predicts that the disk extends down to the last stable orbit at 6R g . From the fits 
to the May (HS/IS) data we see that indeed the reflector subtends a large solid angle with respect 
to the hard X-ray source and the spectrum is strongly relativistically smeared, although the data 
prefer an inner radius of ~ 20 i? g , and are formally inconsistent with 6 R g . Relativistic corrections to 
the illumination emissivity will be important at these small radii, although the results of Reynolds & 
Begelman (1997) suggest that this should lead to a steeper dependence, rather than the flatter one 
needed for the data to be consistent with 6 R g . A more likely possibility to retrieve a disk down to 
the last stable orbit is that the strong ionization of the reflector seen in this state is not adequately 
modeled (the ionization balance derived from Done et al. 1992 is rather approximate), and that this 
has distorted the fit. This strong ionization of the upper layers of the disk probably results from 
the intense gravitational energy release within the SS flow, and it may lead to Comptonization of 
the emerging soft photons, producing the observed spectral shape of the thermal emission (Ross & 
Fabian 1993; Gierlihski et al. 1997b). 

However, the next stage of evolution in the IS does not seem to be properly described by the 
EMN model. Since the hard X-ray flux was roughly constant between May and June the observed 
dramatic decrease of the ionization, from £ ~ 10 4 to ~ 10, must be related to the disappearance 
of the soft component. Consequently, the power law becomes harder, since the availability of seed 
photons for Compton cooling is dramatically reduced. However, the dramatic spectral change was 
not accompanied by any significant change of i?j n , although the best fit value of / decreases by 
~ 0.2 — 0.3. This could be produced by a small increase in Ri n , from ~ 10 R g to ~ 20 R g , depending 
on the emissivity, which is within the error contours of Figure pi It is clearly incompatible with 
the prediction of EMN that i?j n should be ~ 200 R g for the June data since this would produce no 
significant relativistic smearing of the iron spectral features and too small amplitude of reflection. 

Similarly, the properties of the reprocessed component in the July data (LS) are incompatible 
with the EMN computations. While the effect of smearing constrains R{ n to be > 20i? g , the 
amplitude of reflection, / ~ 0.25, requires Ri n to be ~ 30 R g for the EMN geometry (obtained 
using simple Monte Carlo code for photons propagation without further interactions after their 
generation). With Ri n ks 10 R g , as postulated by EMN, the amplitude / would be negligible. The 
values of / inferred by us are also in agreement with those for Cyg X-l and GX 339-4 in low states, 
from both simple fits to the data and detailed modeling of the continuum emission (Ueda, Ebisawa 
& Done 1994; Ebisawa et al. 1996; Gierlinski et al. 1997a; Dove et al. 1997; Poutanen et al. 1997). 
The amplitude / ~ 0.2 — 0.3 seems thus to be a typical value for BHC in the low state, and shows 
that concept of accretion occurring in a optically thin flow from very large radii (~ 10 4 R g ) cannot 
be sustained. 

Thus, the observed transition from high (soft) to low (hard) state does seem to involve a 
retreat of the optically thick material, as in the EMN model. However, the inner disk radius is 



- 7- 



much smaller than the values postulated by EMN in their intermediate and low states. Since their 
treatment of the transition radius marking the edge of the advective flow and the standard SS 
disk is assumed rather than calculated then this may not be an insurmountable problem. The 
major spectral change in their model occurs as i?; n changes from 6 to ~ 30i? g , which corresponds 
well to the observed constraints during the high/low state transition, suggesting that it is indeed 
associated with the innermost regions of the accretion flow becoming optically thin. However, the 
sudden switching of all the accretion flow into an optically thin state over a very small range in m 
is a major facet of their model, and is clearly inconsistent with the observed behavior of BHC in 
general and Nova Muscae in particular. 

This research made use of data obtained from the Leicester Database and Archive Service at 
the Department of Physics and Astronomy, Leicester University. CD. acknowledges support from 
a PPARC Advanced Fellowship. 
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Fig. 1. — Best fit spectra for May 18th data (a) and June 13th data (b) and the fit residuals. 
The upper solid histograms are the total spectra, the dotted histograms show the soft thermal 
components and the dashed histograms show the primary power law spectrum The shape of the 
reprocessed component (lower solid histograms for the reflected continuum and the Fe line), deter- 
mined mainly by ionization parameter £, changed dramatically but its relative amplitude, /, and 
amount of relativistic smearing changed only slightly (see Table [l]) indicating only moderate change 
of the inner disk radius. 
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Fig. 2. — Contours of Ax 2 = 4.61 (90% confidence level for two interesting parameters) as a 
function of the amplitude of the reprocessed component, /, and the inner radius of the disk, i?i n , as 
derived from spectral modelling of the data. Solid contours are for constant ionization parameter 
models while the dashed and dotted contours for the May 18th data are for models with £(r) oc r~ 
and £(r) oc r -4 ' 5 , respectively (see § ||). The results are consistent with a monotonic decrease of / 
and increase of R{ n indicating the retreat of the inner disk as the decline phase was progressing. 



Table 1. Results of model fitting, £(r) = const 



data T a f h C Rin d (R g ) xVdof (spect. bins) 



Jan 11, VHS 


2.03 ±0.16 


30+ 13 

u - ou -0.04 


(1.01S;?) x 10 4 


13+* 


26.1/31 (40) 


May 18, HS/IS 


9 9Q+0.05 
z -^ 3 -0.03 


0.64+.° o f o 


(3.5lf *) x 10 4 




13.5/22 (31) 


June 13, IS 


1.91 ±0.03 


0.57+.° ;f 7 


10+f 


iol 9 4 


20.9/24 (31) 


June 21, IS 


1.83 ±0.03 


0.46 ±0.12 


41" 


13tf 


11.9/24 (31) 


July 23, LS 


1.72 ±0.02 


94+ ' 11 




50t£ 


15.3/24 (31) 


Sept 3, LS 


1.95 ±0.11 


+i.i 


0(f) 




26.9/26 (31) 



a Photon spectral index 
b Covering fraction of the reprocessor 
ionization parameter of the reprocessor 
d Inner radius od the accretion disk 



